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The electromagnetic form factor of the pion is calculated in the "point-form" of relativistic quan- 
tum mechanics using simple, phenomenological wave functions. It is found that the squared charge 
] radius of the pion is predicted one order of magnitude larger than the experimental value and the 

asymptotic behavior expected from QCD cannot be reproduced. The origin of these discrepancies 

■ is analyzed. The present results confirm previous ones obtained from a theoretical model and call 
' for major improvements in the implementation of the "point-form" approach. 

a 

a 

I. INTRODUCTION 

I/-) \ The pion charge form factor has been the object of many studies 0, 0, S 0, 0JE 0J1J3| • The simplicity of the system, 

■ a quark - anti-quark bound state, and the availability of experimental data |lClllll ll2| make it very attractive to test 
both physical ingredients and methods employed in its calculation. In view of a strongly bound system (if one assumes 
constituent quark masses of the order of a few hundred MeV) and the large Q 2 at which the pion form factor has 
been measured, a relativistic calculation is most likely required. There are various approaches to implement relativity, 
ranging from field theory to relativistic quantum mechanics. Among the latter ones, calculations of the pion form 
factor have been done in instant and front forms OlQilEElSIHGiB- Only one has been performed in the less known 

\ "point-form" |l3j . A good agreement with experiment was claimed by the authors Q. 

In this "point-form" approach the parameter determining the elastic form factor is the relative velocity of the initial 
O ,' and final states. Thus, in the Breit frame, the form factor depends on the momentum transfer Q 2 only through the 
quantity, v 2 = Q 2 /(AM 2 + Q 2 ), where M is the total mass of the system. This has the immediate consequence that 
the squared charge radius of the system scales like 1/M 2 , and therefore tends to oo when M goes to zero. Accordingly, 
^ , the charge radius increases when the system becomes stronger bound, contrary to physical intuition. For the pion, 
this argument suggests that the squared charge radius should be of the order of 3/(2 M 2 ) ~ 3 fm 2 , which exceeds the 
experimental value by almost one order of magnitude. In contrast to this result, the calculation of Ref. 9] shows nice 
agreement. Furthermore, these authors show that, in the case of an infinitely bound system, the form factor could 
scale like 1/Q 2 for large momentum transfers, in agreement with the QCD expectation (14L fl5| - However, the general 
derivation of this result does not imply any assumption on the binding of the system. 

Motivated by the above observations, we re-examine in this letter the "point-form" calculation of the pion form 
factor. While doing so, we also have in mind the failure of the approach when applied to a system made of scalar 
constituents [Til Ir7| . The main goal of this study is a clarification of the above-mentioned peculiarities connected 
with recent implementations of the "point-form" of relativistic Hamiltonian mechanics 1 . A deeper understanding of 
this relatively unknown approach should certify that the method is on a safe ground. With this aim we will use simple 
wave functions (a Gaussian one as used in Ref. [!j and a Hulthen one), which allow analytic calculations, providing a 
better insight on some of the features that the approach evidences. 
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1 As it has been observed in our recent work lift Il7l 1 1 Hi and in Ref. Il9l , this implementation is not identical to the form proposed by 
Dirac L'1.1 in that it does not involve a quantization performed on a hyperboloid. To emphasize this difference, we will put the expression 
"point-form" between quotation marks throughout this paper. 
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II. EXPRESSIONS OF FORM FACTORS IN "POINT-FORM" 

The expression of the pion form factor can be defined quite generally as: 

y f2ETlE l (f\J»\i) = F 1 (q 2 ) (P; + If). 



(1) 



By identifying the matrix element of the current with its expression in terms of wave functions, an expression for the 
form factor Fi(q 2 ) can be obtained. We first remind below results obtained for spin-0 constituents and then present 
expressions for the spin- 1/2 case. 

A. Spin-0 constituent case 

The matrix element relative to the "point-form" form factor for a system of spin-less particles may be written: 

1 



yj2E f 2Ei {f\J"\i) = y/2M f 2M 4 ^3 J d 4 P d 4 Pf d 4 Pl d Vf d Vl 9{X f ■ Pf ) 6(X f ■ p) 9(X Z • p) 9{\ ■ Pl ) 
xS(p 2 m 2 ) S(p 2 - m 2 ) S(p 2 - m 2 ) 5\ Pf +p A/7?/) S 4 ( Pl +p X m ) <j> f ((^r^) 2 ) <k ((^) 2 



^(p f +p) 2 (pr+p) 2 (P^f+Pt) 



(2) 



Notations have been explained in Ref. ^(| . Let us mention here that the four- vectors A represent the velocities of the 
system in the initial and final states, X^f = P^/M. 




FIG. 1: Graphical representation of a virtual photon absorption on a pion together with kinematical definitions. 

The kinematics is otherwise given in Fig. ^ The 77 variables have been introduced to make the covariance explicit. 
After performing some of the integrals, the matrix element of the current may be written: 

i_ r dp h ((^) 2 ) 4* ((^) 2 ) 

(2tt) 3 7 2e p (2X i -p i ) (2X r p s ) 



^2Ej2E~ {f\J»\i) = v /2M / 2M i \ ^ 



\J(p f +p) 2 ( Pl +p) 2 



(3) 



In this expression, it is understood that the arguments p^ and p^ must be replaced by their value that stem from the 
5(- ■ •) functions appearing in Eq. (J2J, which after integrating over the r\ variables imply: 



p. 

Pi +P = (ei + e p ), 

-CJi 



p f +p= — (e/ + e p ) 



(4) 



These relations characterize boost transformations in the "point- form" implementation proposed in Ref. |l3j |. They 
differ from those used in instant- or front-form approaches. 

In the Breit frame, the above expressions allow one to write the form factor, F\(q 2 ), as: 



Fi{q 2 ) 



1 + v 2 



v / T^2 J (2tt) 



dp 



<f>f(Ptf) &(p«), 



(5) 
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where Q 2 = -q 2 , v 2 = Q 2 /(4M 2 + Q 2 ) and p 2 lf = p 2 x +p 2 + (p z ±ve p ) 2 /(l - v 2 ). It can be checked that this last 
expression of the form factor, or Eq. after an appropriate change of variable, leads to the relation: 

f iik 

Fi( ° )= y j2^f 4>2{k)=h (6) 

in complete agreement with the charge of the system on the one hand, the ortho-normalization of the solutions of a 
mass operator on the other hand (no extra k-dependent factor is introduced in the integrand). 

The expression of Eq. may not be unique however. Due to the somewhat arbitrary appearance of the mass 
term, M, at the r.h.s., which should come out dynamically an improved expression was proposed in Ref. Apart 
from the fact that it does not allow one to get analytic formulas, it does not add anything to the present discussion. 



B. Spin-1/2 constituent case 



In order to calculate the pion form factor, one has to take into account the spinorial properties of the constituent 
quarks as well as the pseudo-scalar nature of the pion. A minimal covariant factor accounting for these features 
involves the following trace of 7 matrices: 



J" 



4Tw 75 (m+ :-^ ^ 



(m + j-pf) (m — 'j-p) 

n 75 ^ 



Pi (Pf -P) + Pf (Pi ■P)-P l " {Pr -Pf)+m 2 (pf + p>} + p" ) 



(7) 



At zero momentum transfer and when all particles are put on mass shell, this quantity becomes identical to 
y/(pf + p) 2 {Pi + p) 2 (p 1 } + Pi), which is just the factor that appears in the expression of the matrix element for 
scalar constituents, see Eq. l|2"jl. To get an expression for the pion form factor consistent with the normalization 
adopted for the solutions of the mass operator, it is therefore sufficient to replace this factor for scalar particles by 
the above expression / M . One thus gets the following matrix element for the J^ 1 current: 



y/2E f 2Ei (f\J>*\i) = yjTMf 2M ( 



(2tt) s 



d A pd A pf d 4 pi drjf dm 0(A/ • pf) 6>(A/ • p) 9(\i • p) 9(\i ■ pi) 



x 8{p 2 



S(p 2 f 



5(pl 



2 )S 4 (p f+P - Xfttf) 5\p l+P - <p f ( 



Pf-P-,2 



P f 



x 2 



Pt (Pf ■ P) + Pf (Pi ■ P) - P 1 * (Pi ■ Pf) + m (pf +Pf+ P 11 ) ■ 
In the Brcit frame, the matrix element of Eq. (JSJ leads to the expression for the form factor: 

dp 



Fi(Q 2 ) = Vl^ 



(2ny 



<t>f{Ptf) <t>i(pu). 



(8) 



(9) 



It can be checked again that this last expression of the form factor is in complete agreement with the charge of the 
system and the ortho-normalization of the solutions of a mass operator. 

It is noticed that the form factor of Eq. 10 differs from the one obtained with scalar constituent particles by a 
factor (1 — v 2 )/(l + v 2 ), see Eq. |J5J. This last factor is close to 1 — Q 2 /(8m 2 ) for small momentum transfers and 
small binding energies (M ~ 2 m). It may be identified with the Darwin- Foldy term sometimes introduced in the 
calculation of form factors involving spin-1/2 constituents. A similar factor can be found from Ref. |2j as a result of 
Melosh rotations where, however, a quantity 2m appears in place of the pion mass in the definition of v. 



III. ANALYTIC RESULTS FOR FORM FACTORS 



In a few cases, analytic expressions of form factors given by Eq. @ can be obtained. This is the case for a Gaussian 
wave function, that has been used widely in the field (Refs. an< ^ a Hrdthen one 2 . The latter one is of 



2 The relevance of power-law wave functions with respect to hadron form factors was especially noted in Ref. 
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interest because it leads to the correct asymptotic power law for the form factor in the case of scalar constituents. 
These wave functions with the appropriate normalizations are given by: 



^V ^ |^exp[-fc 2 /(2 6 2 )], 



<j> H (k) 



'An 



^2a0(a + 0f 
(a 2 + k 2 ) (0 2 +k 2 )' 



(10) 



The parameters may be chosen as free parameters or fitted to the matter radius. For the Hulthen wave function, a 
is in principle related to the binding energy by a 2 = m 2 — M 2 /A, only is really a free parameter. With these wave 
functions the integral in Eq. may be calculated analytically. One gets in the Gaussian case: 



Ff(Q 2 ) = Vl-v 2 / dp 



1 



exp 



I f 2 pfi+v 2 e 2 



p 



y/2 ^3/4 

and the associated charge radius and asymptotic behavior are given by: 
dFf (Q 2 ) 



-6 



1 + ^-1. *W 



(1 - v 2 f 

VTT^ 2 



1 16 Mt 



exp 



dQ 2 AM 2 

In the case of the Hulthen wave function, the form factor reads 

dp 



Q2 ^°° s/2 Q 4 



exp 



l-v 2 b 2 



Q 2 Am 2 



16b 2 Ml 



Ff (g 2 ) = vr^ 



8ira0 {a + 0)" 



(11) 



(12) 



y/ m 2 — 2 



\J m 2 — f3 2 



The charge radius and high Q 2 properties associated to this form factor are given by: 



(r 2 ) H = 



6 



a (70 2 + A0a + a 2 ) m 2 - a 2 4 + 5 3 a + 12 2 a 2 + 5 a 3 + a 4 



AM 2 V 3/3 2 ((3 + a) 

P^oo 3 1 

~^ AG 8c?' 



2 a 2 (0 + af 



F»(Q 2 ) 



Q 2 ^oo 



64 Mil a (0 + a) 2 p^oo 16 Mi arctan[Vm 2 - a 2 /a] 

a - 



Q 8 (to 4 + . . .) 



Q 



V to 2 — a 2 



(« 2 +?i + (^) 2 ) (/? 2 +pi + (^) 2 )(« 2 +pi + (^) 2 ) (p 2 +pl + ( p -^y 

(1 — v 2 ) 2 a0(a + 0) ( arctan[t)VTO 2 — a 2 /a] arctan[2 v\Jw? — a 2 /(a + (3 — v 2 (f3 — a))) 
v (a-0) 2 y Vto 2 - a 2 \Jm 2 - a 2 

arctan[w \/m 2 — f3 2 / 0] arctan[2 v\/m 2 — 2 /(a + + v 2 (0 — a))] \ 



(13) 



(14) 



where the dots in the last expression represent contributions that can be neglected for finite values of 0. An interesting 
limit is the Coulombian case (0 = a), considered in Ref. ^(| for scalar constituents. The expressions in Eq. (|14f) then 
simplify to read: 



(r 2 Y 



6 

M 2 



256 Mg a 4 
Q s TO 4 



4 a 2 



Ff(Q 2 ) 



(15) 



IV. DISCUSSION 
A. Charge radius 

As it can be observed from Eqs. |Q and ijTill . the squared charge radius scales like the inverse squared pion mass, 
confirming the hand-waving argument given in the introduction and, at the same time, the consequence that the 
charge radius increases with increasing binding energy! 

The second observation is that the squared charge radius has a minimum value, as a function of the open parameters, 
that depends on the pion mass. It slightly differs with the model (15/(4 M 2 ) for the Gaussian wave function and 
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3/Af 2 for the Hulthen one). However, even this minimal values, (7.5 fm 2 and 6fm 2 ) exceed by far the experimental 
one, 0.43 fm 2 [llj. Independently, the contribution of the bare matter radius of the pion to these values is not small. 
This part, 3/(86 2 ), is enhanced by {2m/M r ) 2 , which amounts to a factor 8-18 for quark masses in the range of 0.2 
- 0.3 GeV. For a squared matter radius of 0.05 fm 2 , of the order of what is generally predicted by quark models, the 
corresponding squared charge radius is already of the order of the experimental one. The existence of a minimum 
value largely exceeding the experimental one has the striking consequence that there is no way to fit the parameters 
of the model to reproduce the experimental pion form factor. 



B. Asymptotic behavior 

The asymptotic form factor for a system made of spin-1/2 constituents is expected to scale like Q~ 2 0, In 
addition, the coefficient of the Q~ 2 factor should be proportional to the strength of the short range quark-quark 
interaction (up to log terms) and to the square of the wave function at the origin. These expectations are independent 
of any assumption on the mass of the constituents or on the size of the system. 

Evidently, the Gaussian model cannot reproduce the correct power law behavior of form factors in general. Only 
in the zero-size limit, when b 2 — > oo, one obtains iq cx (M 2 /Q 2 ) 2 . Apart from the fact that the power of Q is not 
the right one, the overall coefficient does not scale like the strength of the quark-quark interaction. This illustrates 
the statement that reproducing the QCD power-law behavior with Gaussian wave functions and standard currents is 
most surely the consequence of a mistake in the theoretical approach employed to get the result. 

Not surprisingly, form factors calculated with Hulthen wave functions provide a power law behavior (Q -8 ), which, 
however, differs from what is expected. For a large part, this is in relation with the observation made in Ref. 21] that 
the transferred momentum to the struck particle in the "point-form" approach is Q 2 [1 + Q 2 /(4M 2 )] rather than Q 2 . 
A better result is obtained in the limit (3 — > oo, in which case the wave function, Eq. I|10|) . scales like k~ 2 instead of 
fc~ 4 . This limit however corresponds to a zero-range force, which is difficult to imagine for an effect dominated by 
one-gluon exchange. 

Apart from that, one can check that the expression of the form factor, Eq. (|13fl . has at least some of the appropriate 
properties which can be easily checked in the Coulombian case, Eq. 115|) . The coefficient in the expression for the 
asymptotic form factor, a 4 , splits into a factor (a 3 / 2 ) 2 , representing the square of the wave function at the origin, 
and a factor a which involves the strength of the short-range force (a = ma s /2). The same observation applies to 
the Hulthen wave function. 

Concerning the asymptotic behavior itself, one generally refers to the Q~ 2 power law from QCD but it seems 
that this is not what one should expect from a relativistic quantum mechanics calculation based on a single-particle 
current. Calculations in the light front approach [2^, in the instant form and in the present work (after disentangling 
the effect of the modified momentum transfer discussed above), rather indicate a Q~ A asymptotic behavior, like for 
scalar constituents. One has to rely then on two-body currents to get the correct Q~ 2 behavior. It can be checked 
that this behavior is obtained without any problem from a single particle current in the case of a scalar probe. This 
indicates that the different power law obtained for spin-1/2 constituents has its origin in the specific form of the 
coupling between quarks and photons (see below). 



C. Discrepancy with an earlier calculation 

Results presented above significantly differ from those obtained in Ref. [9j, where the charge radius was correctly 
reproduced. For an unknown reason 3 , a factor M I /(2m) was introduced in the relation between the momentum of 
the pion in the Breit frame and the momentum transfer, Pij = ±(Q/2) x M,/(2m) (see Eq. (14) of Ref. 9]). The 
velocity v 2 , which is the relevant parameter in the present "point-form" approach, thus becomes: 

2 — Q 2 2 _ Q 2 

~ 4M2 + Q2 U ~16m 2 + Q 2 ' 1 J 

With the last expression, all problems related to the limit M v — > vanish. There is no more lower bound for the 
squared charge radius that prevents one from making a fit of the form factor to experimental values. While the above 



3 This factor was claimed to be necessary to get the non-relativistic limit right but the sense of this statement is not transparent and the 
recipe was not used anymore in later works. 
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modification brings the pion form factor closer to those obtained in other forms, there is no justification for this 
change within the "point-form" approach. Evidently, the more reasonable character of the results obtained with this 
correction suggests that there may be some truth in it [l8j but no tentative explanation was presented by the authors. 

The other difference concerns the asymptotic form factor obtained with a Gaussian wave function. We gave 
arguments to discard the use of such wave functions for looking at the asymptotic behavior of form factors. It 
however remains to explain why, in the zero-size limit, a Q~ 2 behavior was obtained in Ref. |9j while we get Q~ . 
This is most likely due to the matrix element of the current, which in the Breit frame reads: 



(m + e Q/2 ) 2 ~ m + e Q /2 



U 7° U ~ 1 ^ , J ; 2 ' ~ . (17) 



Our expression for the form factor implicitly accounts for the full structure of the spinors, thus evidencing some 
dependence on the momentum transfer, as can be seen from Eq. I|17|l . This tends to reduce the value of the matrix 
element of the charge density. In Ref. this matrix element was assumed to be a constant: a[ \ J° (0) \p± a\) = 
e i &o' cri j which differs from Eq. (|17fl by a factor (m + e.Q/2)/ (2 m). It can be checked that a full account of Eq. i|17fl . 
together with the redefinition of v 2 as given in Eq. i|16|) and the modified relation between the momentum transfer 
and the momentum of the struck particle in the "point-form" approach, provides an extra factor (1 — v 2 ). Dividing 
Eq. (jl 1|> by this quantity, the form factor would read: 

Fl (Q 2 ) = ^=£=, (18) 
V 1 + u 

which is close to the results of Ref. [jj with a quark mass fitted to the experimental charge radius (to = 0.22 GeV). 



D. Numerical estimates of form factors 



In Fig. [21 we provide some numerical estimates for form factors calculated with Gaussian and Hulthen wave functions. 
This is done for what we consider the most preferable case with respect to a reasonable choice of parameters and 
a favorable comparison to experiment. For the Gaussian wave function, the parameter b is chosen such as to get a 
matter radius of 0.025 fm 2 , which is at the lower bound of predictions [23, giving b 2 = 0.6 GeV 2 . The quark mass is 
taken as m = 0.22 GeV. For the Hulthen wave function, we take the limit (3 — > 00 and the parameter a is fixed to 
give the same matter radius as in the Gaussian case, giving a = 0.2 GeV 2 . The quark mass here is chosen such that 
m 2 -a 2 = M 2 /4. 

In addition, Fig. [21 contains the form factor given by Eq. I|18[) . with the redefinition of v 2 as discussed before, see 
Eq. (|16|) . With a value of m = 0.22 GeV, one obtains a good representation of both experimental data and results 
obtained in Ref. 9] for similar conditions. Evidently, the Gaussian and Hulthen form factors miss the experimental 
values by a large factor. In the very extreme but unphysical limit of a point-like pion, the following upper limits are 
obtained: 

F?m> = ^=3f , W)> = log^- (19) 

VI + v 2 2v 1-v 

These ones, also reported in Fig. [21 are still far below the experimental values. 



V. CONCLUSION 

In a previous work [l6l Il7| , "point-form" form factors were calculated and found to evidence a large discrepancy, 
both at low and high Q 2 , with what could be considered an "exact" calculation. The system under consideration was 
academic however. In the present work, we applied the same approach to the pion form factor whose low Q 2 behavior 
(related to the charge radius) and high Q 2 behavior are determined sufficiently well from experiment to make reliable 
statements. Results we obtained, based on a single-particle current, evidence quite similar features: too fast fall off 
at both low and high Q 2 . In the first case, they point to a squared charge radius more than one order of magnitude 
larger than the measured one. In the second case, the asymptotic behavior misses by many powers of Q the l/Q 2 
power law expected from QCD. 

For our purpose, we used simple wave functions which offer the advantage of providing analytic results. Better 
wave functions could be used. However, the comparison with results from Refs. |l6l Il7| does not show qualitative 
differences as far as the most striking features are concerned. At low Q 2 , results obtained with Gaussian and Hulthen 
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FIG. 2: Pion form factor calculated with different wave functions. The left hand graph shows the low-Q 2 behavior in linear 
scale, the graph on the right hand side shows the form factors, multiplied by Q 2 , in logarithmic scale up to Q 2 = 10 GeV 2 . The 
different curves correspond to Eq. 1111 (solid), Eq. 1131 with 8 — » oo (dotted), Eq. I19H (long- and dot-dashed) and Eq. 11811 
(dashed), respectively. Experimental data are from Refs. [TH. ITU ll2T|. 

wave functions are very similar. Both point to an uncompressible value of the squared charge radius of the order 
of 6fm 2 , an order of magnitude larger than the experimental value. This prevents one from making any fit to the 
measured form factor. At high Q 2 , the model dependence may be more important but, apart from the fact that the 
correct power law behavior is missed, the most favorable cases, such as the zero-size limit, are physically irrelevant. 
In all cases we considered, the "point-form" form factors are far below experimental data. 

Present results differ strikingly from those obtained in Ref. For the largest part, this is due to a different relation 
between the pion Breit-frame momentum and the momentum transfer which was used by these authors, allowing them 
to escape the constraint of a minimum value of the squared charge radius determined by the pion mass. However, 
this relation has no theoretical support within the "point-form" approach. 

On the basis of a theoretical model, it was concluded in a previous work that the present implementation of the 
"point-form" approach together with a single-particle current requires major improvements |l6j . The present results 
for a physical quantity, the pion form factor, lead to the same conclusion. Large contributions of two-body currents 
are needed [24[ . An alternative would be to improve the present implementation of the "point-form" approach |2^| . 
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